Dynamics of Diblock Copolymers in Dilute Solutions 



Radu P. Mondescu and M. Muthukumar 

Department of Physics and Astronomy, and Polymer Science & Engineering Department and Materials Research Science and 

Engineering Center 
University of Massachusetts, Amherst, MA 01003 

In the present work we consider the dynamics of freely translating and rotating diblock (A- 
B), Gaussias copolymers, in dilute solutions. Using the previously developed multiple scattering 
techniqueou — applied to the study of suspensions of spheres and polymers — we have computed the 
diffusion and the friction coefficients Dab and Cab, respectively, and the change <5?7ab in the viscosity 
of the solution as functions of x = and t — jj, where Na, N are the number of segments of the 
A block and of the whole copolymer, respectively, and Ia, Ib are the Kuhn lengths of the A and B 
blocks. Specific regimes that maximize the difference between the diffusion constants of copolymers 
with distinct t values, which lead to increasing the efficiency of separation processes, have been 
identified. 



I. INTRODUCTION 



Currently, due to both the practical significance and the theoretical challenge, there is a sustained interest in 
the study of transport properties of copolymers in solution, mainly in relation to electrophoresis and associated 
phenomena. Basically, one wishes to tag the macromolecules of interest with a neutral chain to make the friction 
coefficient of the resulting object molecular weight dependent for all measurement conditions, in order to increase the 
efficiency of the separation process.Q As there are experimental regimes for which the molecular weight dependence-js 
lost and the separation is very weak if not impossible (for a recent calculation for polyelectrolytes see MuthukumarO) , 
a broader understanding of the dynamics of composite objects in solution becomes critical. 

To gain insight in the dynamical behavior of an arbitrary heterogeneous polymer in solution, we will address a 
simpler, nevertheless illuminating problem: consider two linear, Gaussian polymers A and B, with lengths/number of 
segments La/Na and Lb/Nb, respectively, and Kuhn lengths Ia and Ib, joined to form a composite chain (A-B) of 
total length L = La + Lb and total number of monomers N = Na + Nb ■ In this paper we calculate the diffusion 
and the friction coefficients of the diblock copolymer A-B and the change in the viscosity of the solution due to the 
copolymer in dilute solutions, as functions of the dimensionless variables t = and x — 

Our calculations are based on the cluster expansion theory and on the effective medium hypothesisJ3'S~i The main 
simplifying assumption is the absence of excluded volume interaction. Other approximations and features of our 
model are: 

a) no interaction between the polymer chains except the hydrodynamic interaction is present. 

b) the solvent is incompressible and the steady state limit is considered. 

c) the solution is assumed infinitely diluted (one chain limit). 

d) the usual preaveraging approximation,^ where one replaces the Oseen tensor G(r,r') with its configurational 
average (G(r, r')), is employed. 

e) the hydrodynamic interaction is enforced by using stick boundary condition. 

The model is not limited by the Gaussian nature of the two polymer components, because the Kuhn lengths Ia and 
Ib could be effective quantities modeling chains with excluded volume and Coulombic interactions (assertion valid in 
the hydrodynamic limit; at short distances the fractal structure of the polymer becomes transparent). Also, even if 
the diblock- copolymer could be regarded as a particular example of a branched polymer, extensively studied in the 
literature )13't2l the important distinction in our work is that we allow for unequal Kuhn lengths for the two branches, 
with interesting consequences in the calculations — the structure factor of the copolymer must be computed now using 
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a bivariate Gaussian probability distribution function. Finally, the present approach allows for obtaining exact results, 
albeit numerical except some limits, and no approximation of Kirkwood-Riseman type is necessary. 

Concerning the outline of the paper, we next develop schematically — emphasizing the new elements — the underlying 
theory, then we show explicitly the calculation of -Dab, Cab and of <5?7ab- We conclude with a discussion of the results. 



II. CALCULATION OF THE STEADY-STATE VELOCITY FIELD: INFINITE DILUTION LIMIT 



The dynamical properties of the system can be deduced from the knowledge of. the average velocity field at any 
point in the solution. As this type of computation is already very well documented^ we will show detailed derivations 
only when necessary. 

Briefly, the main steps of our computations are: first we write the velocity field in the solution in the presence of the 
copolymer; we apply the boundary condition to eliminate the unknown friction forces between the polymer and the 
fluid and then we average over the random position and the configuration of the chain; next we construct the cluster 
expansion in terms of the self-energy (Hartree-Fock diagrams) £ of the fluid and we express it as a function of the 
known (measurable) physical parameters of the system; finally, we use the self-energy to extract — both analytically 
and numerically — the quantities of interest. 

Assuming that the diblock copolymer A-B with total number of segments N is immersed in an ideal, incompressible 
fluid obeying the linearized, time independent Navier-Stokes equation (with the solvent density equal to one), we can 
write the equation followed by the velocity field v(r) everywhere in the suspension as 

- 77oAv(r) + Vp(r) = F(r) + / ds 5(r - H s )a s , (2.1) 

Jo 

where F(r) and p(r) are some external force producing the flow and the pressure, respectively, ijo is the kinematic 
shear viscosity coefficient and <x s = er(s) is the force exerted by the bead located at R s = R(s) upon the fluid, s being 
the arclength al ong the chain (for notational convenience we will subsequently write s as a subscript). The formal 
solution of Eq. ( |2.l[) is readily given by the integral representation (equivalent to the single layer solution of standard 
boundary methods) 



v(r)=v (r)+ / ds G{r-R s ) • <x s , (2.2) 
Jo 

in which vo(r) = fdr' G(r — r') ■ F(r') is the externally imposed velocity field and G is the Oseen tensor: 

1-kk 



G(k) = 



T] k 
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Here, 1 is the unit tensor and k is the unit vector pointing in the k direction and the Fourier transform is given by 

/d\i 
-_G(k)expHk.r). (2.4) 

As mentioned, the dynamics of the copolymer chain and of the fluid are coupled through the no-slip boundary 
condition 

R s =R + u; x S s =v(R s ), (2.5) 

where the velocity of the center-of-mass Ro is constant and we further assume that the angular velocity u> is constant 
and equal to the configuration averaged angular velocity of the chain. To this equation we attach the physical 
constraints of force-free and torque-free motion (neglecting the inertial terms) of the diblock chain: 
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— f (total force) = / ds cr s — 
Jo 

— M (total torque) = / ds S s x er s = 0. (2.6) 
Jo 

Above we used the relation R s = Ro + S s , where S s represents the position vector of the monomer at s about the 
center-of-mass (Ro) of the copolymer. 

Next one needs to eliminate the unknown forces <x s , a task accomplished by using the boundary condition ( |2.5| ) 
and defining the single-chain inverse operator G _1 as 



/ ds' G 1 (R s , R fl , ) G (R y - R s ,, ) = S ss „ 1 . 
Jo 



We obtain then: 



(2.7) 



a s = dfl / G- 1 (R a ,R a /)-(Ho + wxS^-v (R^)), (2.8) 
Jo 

where the variables s, s' should be read as the arclength arguments of R(s), R(s'). Inserting er back in ( |2.6| ) we can 
calculate the total force f and torque M acting upon the polymer as 

- f = g T ■ Ro + Str ■ w - J J ds ds' G-) • vo(s') (2.9a) 

-M = gRT ■ R + g R ■ u> - J J ds ds' S s x G-) • vo(s'), (2.9b) 

where we manipulated G~, as a dyadic tensor. The newly introduced tensors gr.R provide, upon averaging over the 
distribution of the segments s and s' , the translational and rotational friction coefficients of the copolymer: 

Cab = (gr) = / JJ^ dsds ' G ss^J (2.10a) 

Cab = (gfl) = -(^JJ ds ds ' Ss x G -' x S *' y ( 2 - 10b ) 

The remaining coefficients Qtr.rt (the cross translational-rotational and rotational-translational terms) vanish upon 
preaveraging, so will be further discarded. For future reference, their expressions are: 

L 



Etr = ~ J J dsds' G' 1 , x Sy (2.11a) 
Ert = II dsds'S s x G-). (2.11b) 



Returning to the constraints (2.6), we now eliminate the unknowns Rq and u> using 



Ro = g^ 1 • J J ds ds' G~) ■ v (s') 
w = gfl 1 ' J J dsds'S s x G~ s ) • v ( S '), (2.12) 

where the inverse operators g^ 1 and g^ 1 are defined by the relations &p l R ' St.r. = 1- Inserting the values obtained 



for Rq and u> in Eq. fl2.8|) we find that 
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J\s> Ui-) fl Q L dtdt' • g^ 1 • G-,], 

dtdt' G~ t l x S t ■ g^ 1 ■ St x G^ 1 ,} • vo(s') (2.13) 



By introducing this expression back in Eq. (2.2) we get the solution for the velocity field v(r) 



v(r) = v (r) + / / / dv'dY"dv"' G(r - r') • T(r', r") • v (r"). (2.14) 



T, the flow propagator that embodies the effect of the polymer chain upon the solution, is given by 

i-L 



T(r,r') =~JJ dsds' S(v - R s ) T ss , S(v' - R s ,) 

T ss , = T(R s ,R y ) = G~) -11 dtdt 1 G^ 1 ■ g T 1 • G-,l, (2.15) 



L 



dtdt'G^ 1 xSfgx 1 -S t , xG-,1,, 

where, as previously mentioned, the subscripts s,s',... represent the arguments R S ,R S ',... of the operators. T 
depends in general upon the position, the structure and the geometry of the copolymer. 



Finally, we obtain the equation for the macroscopic velocity field by averaging Eq. (2.14) with respect to the random 
distribution of Ro : 

u(r) = <v(r)) = v (r) + G * (T) * v | r , (2.16) 

where we introduced the notations * and | r to denote the convolution operation and its final argument: A * B| r = 
/dr'A(r,r') • B(r'). 

The connection with the experimentally measurable quantities is made by introducing the self-energy tensor S(r, r') 
defined by 

J dr 1 S(r,r') • (v(r')) = S * (v) = I J dsS(r — R s )cr s \ . (2.17) 



Next we average directly Eq. (2.2) over Ro and we insert the expression of X from above, which yields 

u(r) = v (r) + JJdv' dr" G(r - r') • S(r', r") • u(r"). (2.18) 

Taking into account that vo(r) = G * F| r , we can formally manipulate this equation to obtain 

u(r) = g *F\ r = Jdr'g(r-r')-F(r'), (2.19) 

in which Q — the effective Oseen tensor — is given by 

0(r,r') = (l - G*S)- 1 * G| (rjr/) 

^--^tm- <-») 

Here, by taking the inner-product with the projectors 1 — kk and kk that obey (1 — kk) • (1 — kk) = 1 — kk and 
kk- (1 — kk) = 0, we decomposed the tensors in their transverse and longitudinal parts as in A =Aj_(l — kk) + Ay kk, 
with A a symmetric tensor. After Fourier transforming and projecting the transversal and longitudinal components 
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of the tensors, we ca n calculat e the self-energy as a function of the flow propagator T by equating the expressions of 
u(r) from Eq. ( [2. 16 ) and Eq. ( 2.18 ). These operations yield 



S(k) = 



Tl(*0 



l + G±(k)T x (k) 



(i-kk) + T|,(fc) kk, 



(2.21) 



where T(k) is the Fourier transform of (T(r — r')) and Gy = (Eq. (2.3)). Note that the equation for £ is exact. 
In the dilute limit we work with, we can further approximate 



E(k) ~ T(k), 



(2.22) 



where pre-averaging upon T(k) is understood. 

From the structure of the effective Oseen tensor (2.20) we identify the change <5?7ab in the viscosity of the solution 
due to the added copolymer, in the hydrodynamic limit |k| — > 0, as 



r] r] k^o ok 1 



(2.23) 



III. CALCULATION OF Dab, Cab AND ^ab 



In this section we present the detailed derivation of the transport properties of the system. The essential aspect is 
the computation of the structure factor, which involves, for the diblock copolymer case, the use of a bivariate Gaussian 
probability distribution function. This calculation is presented in Appendix [a|. For any operator A(s, s') depending 
on the arclength variables s,s' we will use the following Fourier representation (Rouse modes expansion in polymer 
language) : 

oo oo 

A(s, s') = jj ]T £ A pp ,e^ s -^ V (3.1) 

p—~ oop' — — OO 

A pp , = JJ^ dsds' A{s lS ')e-^ ps+ ^ p ' s ' , 

where we remind that L is the total length of the polymer chain. Also, the inverse of A is defined by: 

S(s - s") 1 (3.2) 
L 2 <L C » 1 



/ dj , A- 1 («y) • A(s'.s") 
Jo 

oo 

V A" 1 -Air, 



A. Diffusion coefficient Dab 



From Eqs. ( 2. 9a ), ( 2. 10a ) we obtain immediately — after preaveraging (which is equivalent to replacing all instances 
of the Oseen tensor G(s, s') with its configurational average (G(s,s'))) — that Dab, a tensor quantity in general, is 
just the inverse of the friction coefficient Cab (Einstein relation). Thus we write 



° U; D AB 1 = -^ 1 1 'I-''-' (G(.s..s'i) 



k B T 



(3.3) 



L 2 



Go 
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where we used the Fourier expansions (2.4), and the integral J dflk (1 — kk) = Splitting the double integral 



over the arclength variables s and s' in the regions [0, La]', [La,L], inserting the value of the exponential term (the 
structure factor of the diblock copolymer) from App. ^ and integrating over fc, we get 



16 1 1 1 



Dab = 

k B T iTy/ir L 2 770 IaIb 



1 + (Ri + R 



,3/2 



(3.4) 



where Ra 



LaIj 



LbIe 



are the radii of gyration of the two segments A and B, respectively, of the 



6 , R B - y 6 

copolymer. It is readily checked that .Dab reduces to the classical Kirkwood-Riseman (K-R) resul 



espe 

E 



D 



9^^770 R 3 ' 



(3.5) 



with i = ift)Wl2ir 3 N = 6nV2Trr] a Rg and R g = in the limits l A = Ib, Na = or N B = 0. 

The dimensionless variables that control t he d iffusion of the copolymer are made apparent by normalizing the 
expression (3.4) for Dab to the classical value (|3.5[) of D of a polymer with the length Ljy = NIa, where N — Na+Nb- 



Defining t = -fi- and x 

l A 



N A 
N 



yields, upon normalization, 



L»ab i 2 (l - xf/ 2 {\ - t) - a; 3 / 2 (l -t) + [x + i 2 (l 



13/2 



D 



t[x + t(l -x)] 2 



(3.6) 



which is the expression plotted in Fig. [I] as a function of x and t, where t has values both smaller and larger than 1. 
When t = l,a;=lora; = we obtain the meaningful limits 1, 1 and ~, respectively, the last result (when x = 0) 
representing just the diffusion coefficient of a polymer of l engt h Lb — NIb normalized to the corresponding value of 
a chain with length L^ = NIa- When t = 0, the limit of (3.6) is -4=. From Fig. [l]we notice that a slight increase in 



the monomer fraction of the A block in the A-B copolymer when t = ^- < 1 has a much stronger effect on Dab than 
an even large increase of x in the regime with t > 1. 

A more relevant relation — particularly from an experimental point of view — is obtaine d by computing the ratio 
between Dab and the diffusion coefficient Da of the block A given by the standard result (3.5) (with N replaced by 
Na), which immediately gives 



D AB _ r t 2 {l - x) 3/2 (l -t)- .t 3 / 2 (1 -t) + [x + t 2 (l 
D A ~ t[x + t(l-x)] 2 



13/2 



(3.7) 



This expression is plotted as a function of x and t in Fig. g for a sequence of values of t in the range [1/10 — 10], the 
points with t > 1 being marked by circles joined by lines. Now at each fixed value of x (correspondingly Na) one can 
compare the change in the diffusion coefficient of the diblock copolymer with respect to the diffusion constant of one 
of its components, which has direct implications in separation techniques. 

We remark that at small x values and t < 1/4 a maximum in -^^a. occurs. Also, -Dab exhibits a change in 
curvature when the t parameter goes from t < 1 to t > 1, which motivated us in calculating the separation curves 
A(Dab/Da), defined as the difference between two values of the ratio Dab/ Da corresponding to distinct values of 
t, as a function of x. These curves are displayed in Fig. |[ Physically, the behavior of the diffusion coefficient of 
the diblock copolymer opens the way to control the process of separation of macromolecules by tagging the targets 
with other weight-controlled polymers. In particular, one can think of modifying the diffusion coefficient of a given 
polymer chain A by attaching to it another polymer B, with t having a prescribed value. This is illustrated in Fig. ^ 
by the maxima occurring in the A(Dab/Da) function. There is also a range of the x values that offers an increased 
efficiency of the separation. This domain is much more localized for t < 1 than for t > 1, but the absolute value of the 
peak of A (-Dab /Da) is highest when one compares — at the same x ratio — two copolymers with t greater and smaller 
than one. This is reasonable — t small means I a big, so even small variations in x have more significant effects than 
for pairs of copolymers with t comparable to one. 
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B. Friction coefficient (ab 



From Eq. ( [2. 10a ), expanding the inverse Oseen tensor G~! in double Fourier series ( |3.l|) , we find that £ab = Gqq 1 . 
Because the G pp ' tensor, due to the diblock structure of the copolymer with different Kuhn lengths, is not diagonal 
in the p index, the calculation of the friction coefficient is more difficult. In general, one needs to decompose the G~ p } 
matrix in its diagonal and non-diagonal parts and then to use a Born expansion up to the desired orderu. Fortunately, 
next we will show that the Gqq 1 factor can be computed by direct numerical inversion. 



We start with the expression of G pp ' elements given by Eq. (3.1) 

dk 1-kk ' rh 



G, 



(2tt) 3 T] k 2 



2i-K 



ds ds' (exp[zk • (S s — exp[ ps 

L 



2l7T 



' '1 

p s \ 



(3.8) 



To evaluate the integrals, first we split the domain of integration (0, L) x (0, L) over s and s' in four regions: (0, La) x 
(0, La), (0, La) x(La, L), (La,L) x (0, La) and (La,L) x (La, L), and then we insert the expression of the structure 
factor from App. |Aj. As the remaining integrations are rather intricate, we present the computation of G pp i elements 
in Appendix M. The final result is: 



'pp> \p=£0 



1 



1 



7TV67T770 \JL\P\ 
1 1 



L A 1 
1 



27r 2 V6^?o vxbl *(p'-p) 



La 
L 



1 



a, 



2wiL / 

1 — e ~ 



a ( P ,p'eZ). 



(3.9a) 
(3.9b) 



From Eqs. (3.3;3.4) the Goo element is just L 2 Dab/ '(ksT). The G pp > matrix is hermitian (G pp > = G* p , p ) and 
diagonally dominant, thus securing a well behaved inverse. We remark that only when Ia — Ib, La — or La = L 
the non-diagonal elements of G pp i vanish. 

To further proceed with the numerical inversion, we need to reduce the G pp i tensor to a dimensionless form. This 
task is accomplished by noticing that the Rouse coefficients G pp i can be written as 



G, 



L 



1 



1 



47r 2 v67Tf?o VLIa 



Kpp' (x, t) , 



(3.10) 



where K is a purely numerical tensor, function of the dimensionless parameters x = and t = and with the 
matrix representation (hermitian) given by 



N 



, . 32tt t 2 (l-x) 3 / 2 (l-t)-x 3 / 2 (l-t) + [x + t 2 (l-x)] 3 / 2 
K ° 0{X > t) = ~3 t[x + t{l-x)]>'* 1 



Kpp(£, *)|p^0 

K pp / (x,t)\ p ^o 



4tt 



\P\ 
2 



y - 



1-2/ 



Vt J 



i; y = 



x + t(l- x) 



- 1)[1 -e 2l7r(p '- p)y ] 1, 



(3.11a) 
(3.11b) 
(3.11c) 



where we employed L = La + Lb = NIa(x + t(l — x)). Note that for our choice of t the limit of t — > for K pp / is not 
immediate. Next, applying the definition ([T^) we construct the inverse G pp \ as 



G PP > • Qp'p" - l 2 8 pp „ a 

p' — — CO 



(3.12) 



and substituting G pp i from Eq. ( 3.10 ), we find that 

G" 1 , = 4ir 2 V6^vo \fU~ A K" 1 , (x,t), 



(3.13) 



where the K 1 operator is defined by the equation: 
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Kpp 1 ' (x,t) ■ Vip' P " (x,t) — S pp " 1. (3.14) 

p' — — oo 

E vent ually, recalling that Cab = Gqq 1 and after normalizing to the friction coefficient ( = 97Ty / 7F?7oR.g/4 (see 



Eq. (3.5)) of a homogeneous polymer of length L = Nl&, we obtain 



Cab / 327r 



^ -(a;,t) = — y/x + t(l-x) K^(x,t), (3.15) 

where we dropped the unit dyad from the formulas, as inessential. The limiting cases when the tensor K^ 1 , is diagonal 
can be readily checked analytically. We then get: 

327T 47T 

Koo(x,l) = -5-; K pp (x,l) 



3 ' ^ " 

K 00 (0,f) = ^-^; V(a:,l) = 0forp^' (3.16) 



32tt t ^ , . 32tt 
K 00 (l,i) = — ; K w (x,0) = — 



For these limits the inverse of K o is just Kqq 1 = (K o) 1 and the friction coefficient *&P-(x,t) reduces accordingly to 
the physically correct values: 

^Cc,l) = l; ^(0,t)=t; ^(l,t) = l. (3.17) 

The explicit inversion of K o(a;, t) is carried out by first fixing the values of the parameters x and t. Then we iterate 
the dimensions of the K, pp i(x, t) matrix and we numerically invert it for each \p\ = l,p max until K.Q (x,t) converges 
to the desired precision — we have chosen 10 -7 for computational convenience — at a certain p max . 



The final results for calculated from Eq. ( ft.lE ) have been plotted in Figs. |] and || against the fraction x of A 



beads, for values of t smaller and greater than 1. Points where a numerical evaluation was carried out are marked by 
symbols, their shape corresponding to different values of t. 

An illuminating aspect of the result is that neglecting the off-diagonal components of the Rouse tensor G pp > and 
inverting it directly (then Cab is just L 2 (Goo) _1 ) still leads to the qualitatively correct behavior of Cab, as displayed 
by the continuous lines in Figs. gj. 



C. Intrinsic Viscosity [ijab] 



As seen from Eq. ( 2.23Q , the information regarding the change in the viscosity of the solution is contained in the 
self-energy S(k). To find it we proceed by averaging over the random position Ro of the center-of-mass of the polymer 
and Fourier transforming Eq. ( p. 15 ), which yields 



S(k) ~ T(k) = - I Jd(r - r') J dR exp[zk • (r - r')] 



1 

V 



dsds'6(r - Ro S s )T ss ,S(r' - R - S s ,) 



ds ds' T ss i exp[ik • (S s — S s >)], 



where V is the volume of the suspension. Inserting the full form of T ss / from ( [2.15| ) produces 

S(k) = -i[T«(k) - T< 2 )(k) + T< 3 )(k)], 



(3.18) 



(3.19) 
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with the new factors being given by: 

cL 



T«(k) = J J dsds' G~ s ] exp[zk ■ (S s - Sy)] 



T^(k) = 
T< 3 )(k) = 



ds ds' dt dt' G't 1 • 1 • G~ s \ exp[ik • (S. - S a / )] 
dsds'dtdt'G't 1 x S t • g^ 1 • S t / x G^ 1 , exp[ik • (S s - S s /)] 



(3.20a) 
(3.20b) 
(3.20c) 



Here, we recall that G s ^, is defined by Eq. (2.7) and should actually read G 1 (S S — S s '). From (2.f0) we calculate- 



within the pre-averaging approximation — the frictional coefficients gy and g# as 



( ST ) = fj dsds' (g; s ))=g 



-1 
oo 



(sr) = dsds' (G-)) ■ (S a x 1 x S s 



(3.21a) 
(3.21b) 



dsds' (G- s )).(S s .S s ,) = l± [ Yl E G &'Kp» 



+00 +oo 



3L 2 



-— oo p' — — oo 



where we used (si^ = | (S s • S s ') % that implies ((S s • r)(S s ' • r)) = \r 2 (S s • S s ') for any vector r. Note that 

the double summation over p and p' is equivalent to taking the trace of G^ 1 • F. We also developed all quantities in 
the double Fourier series (3.1) and (3.2). The F elements are computed from 



F / — F /I 

pp — pp 



,ci o \ r 2Z7T 2J7T 

(S s • S s /) exp[ — —ps + 



p's'] 1. 



(3.22) 



We remark that F is a hermitian operator. Averaging S s • S s ' over the distribution of the segments s and s' for the 
diblock copolymer is a rather involved operation and the details are presented in Appendix KJ. 



Next we can simplify the calculations by taking directly the k 2 components of T^(k) in ( 3.20 ). Then TW becomes, 
upon preaveraging,Q 



-fc 2 — 



dsds' (G-}>-<(k-S.)(k-S./)> 



dsds' {G; s ))-(S s -S s ,) = -k 2 ( SR ) 



(3.23) 



where gR is given by Eq. ( 3.21b ). 



(2) 

In a similar manner, it follows that T\ 2 can be written as 



.(2) 



dsds'dtdt' (g;, 1 ) • (g^ 1 ) • (G^ 1 ,) (S. ■ s s 



(3.24) 



1 1 



YJp- fc2 (St 1 ) E G op 1 ' F pp' ' G p'o • 



where the Fourier expansion ( |3.1[ ), (3.2) were once again utilized. For a homogeneous polymer (t = 1), becomes 
diagonal and F o = (App. |c]) so Tl? vanishes as required!. 
Finally, the fc 2 part of T^ 3 ) (k) results from 



,(3) 
■ fe 2 



dsds'dtdt' (G't 1 ) ■ (1 x (k-S s St)) 
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■((k-SyS^xl)-^,) (3.25) 
= -~ (gfi 1 ) k 2 (1 - kk) ■ || ^ dsds'dtdt' (G S Y) 

where we pulled out the g^ 1 factor as it is just a number multiplying the unit tensor and we applied the relation 
(1 x a) • (a x 1) = — a 2 (l — aa), valid for any vector a. Developing all factors in double Fourier (Rouse modes) series 
we obtain 

T$ = - ™4 fc2 ~ ^ ■ E G ^ ' *V E G ,V • F «V> ( 3 - 26 ) 



in wh ich all the summation indices run over the (— oo, +oo) interval. Replacing the sums in terms of the (gi?) coefficient 
(Eq. ( 3.21b| )) and noticing that in our approximation (g^ 1 ) • (gii) = 1, we get the simpler form 



.(3) 
"fc 2 



--(g fl ) fc 2 (i-kk). 



(3.27) 



Collecting all the components of T k 2 from Eqs. ( 3.23| ), ( 3.24 ) and ( 3.27 ) and inserting them in Eq. ( 3.19 ), we find 
that 



S fc2 (k)=-- 



, , too +oo 

ifc 2 (ga ) i - \e ( gT i) ± E G^opj i 



(3.28) 



p=— oo p' — — oc 



\k 2 ( gR ) (i - kk) 



where we used the fact that all quantities are scalar multiples of the unit dyad 1 (e.g. (gR) and (g^ 1 ) are now some 
functions to be computed and G~ p \ , F pp i are the matrix elements of the Fourier representation of G _1 and F). The 



change in viscosity follows from Eq. (2.23) by applying the transverse projector 1 — kk to the previous expression: 



Vo Vr] 



\k 2 (gn) ~ \k 2 (gr 1 ) i Y, G o P % P 'G^ 

p,p> 



(i - kk). 



(3.29) 



Introducing the concentration of the copolymer cab = y jva ' Aav Avogadro's number and Mab the molecular 



mass, we calculate the intrinsic viscosity [t]ab] from 



[VABl 



VAB - T]Q _ Svab VN A y 
cabVo Vo M A b 



(3.30) 



Formally, we have the desired expression. Practically, we need to reduce it to a dimensionless form amenable to 
numerical evaluation. To this e nd, w e write G _1 in terms of the K~ x operator from Eq. (3.13) obeying also (3.14), 
with K(x, t) introduced in Eq. (3.1C). Then we express F pp < as 



(3.31) 



where x = t — l -f- and f PP ' is a dimensionless function of x and t (see App. |^). Inserting back in Eq. ( 3.21b ), we 
obtain the rotationaf friction coefficient as 



(gR, 



+ 00 



= ^-'i/^«) 3/2 L , „ - ^, /2 £ Kpp'Mf^&t) 



[x + (1 - x)tf/ 2 



(3.32) 



p,p' — - — OO 



10 



To check the result we quote some particular values of the f pp i(x,t) function: 

/oo(M) = 0; /oo(0,t) = ;/ 00 (M) = 

fpp(x, 1) = ^. 2 p2 > fpP' ( X i^) = 2tt 2 pp' > /pp( X ' 0) 



(3.33) 



7T 2 p 2 



When i = 1 (the homogeneous case), the K pp / matrix is diagonal, K 1 = (K) 1 and using ( 3.1(f ) (gn) reduces to the 
known result 



24 fm\\ 3/2 ™ i 

(gfl )l t=1 = _ % ^_j ^ 

p— 1 



(3.34) 



Now, as (gy 1 ) is just the inverse of the translational friction coefficient (gr) given by ( 3.21a ), we can calculate 

(gr) = D AB 



t 2 (l - xf/ 2 {l -t)- x 3 / 2 (l -t) + [x + t 2 (l - cc)] 3 / 2 
t[x + t{l - x)} 2 



(3.35) 



from which we recover the Kirkwood-Riseman result when t — > 1, as already discussed. 

Going back to the viscosity formula ( |3.29 ), we use Eqs. ( 3.13| ), ( |3.31 ), (3^21 and ( 3.35 ) to express everything in 
dimcnsionless variables and, after normalizing to the non-free draining resultoH 



6rj _ 1 6 fNl 2 A ^ 3/2 



_ , Zeta(3/2), 

where Zeta(3/2) = Y2pLi (^he Niemann function), we finally obtain 



01] 



Zeta(3/2) [x + i(l - x)] 3 / 2 



+00 



E K-i( x ,t)f; p ,{x,t) 

p,p' — — oo 



64tt t 2 (l - x) 3 / 2 (l - t) - x 3 / 2 (l — *) + [» + < 2 (1 - x)] 3 / 2 



t[a + t(l - a;)] 3 / 2 



+00 



x E K^ 1 (a;,t)/ R p,(a!,t)K^ 1 (a;,t) 



(3.36) 



(3.37) 



Once again , we c an ch eck analytically some limits of this expression by using the matrix elements K pp / and f pp i from 
Eqs. (Ell) and (f]||: 



^f(x,l) = l; iMB(l,t) = l; ^B(0,i)=t 3 . 



(3.38) 



To compare directly the increase in viscosity due to the A -B co polymer to the change due to the polymer A only 
(computed from the K-R result), we only need to multiply ( 3.37 ) with -572 , which then yields 



^ ab (mh4^(m) 



Stja 



£.3/2 fiyj 



(3.39) 



The relative intrinsic viscosity of the suspension follows readily from (3.30) 

[t)ab] _ M A b Stjab 



[rj\ M drj 



(x,t), 



(3.40) 
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where M and [rj] are the molecular mass and the intrinsic viscosity of a solution of a polymer with TV = Na + Nb 
beads and Ia Kuhn-length. 



In Fig. o the viscosity calculated from (3.37) is plotted against x = for different values of the parameter t — — 



From a practical point of view, it is reassuring that, checking once again the validity of the Kirkwood-Riseman 
approximation by computing 8rj ab neglecting all the non-diagonal elements of G pp i , we obtain a very good agreement 
with the exact numerical evaluations. This aspect is illustrated in Fig. ^, where the continuous lines represent the 
calculation done using only the diagonal elements G pp only. 



Based on this finding, in Fig. [7] the viscosity calculated from (3.37) but employing the Kirkwood-Riseman approx- 



imation is plotted against x at various t values. Similar to the observations in the diffusion case, we notice that the 
intrinsic viscosity [^abI/^a] attains a minimum matching the maximum of D AB / D A from Fig. |. This should be 
useful for some technological processes where a control of the viscosity of a solution is desired. 



IV. CONCLUSIONS 



In this paper we have considered the stationary dynamics of an infinitely diluted solution of freely translating and 
rotating Gaussian diblock copolymers. The copolymer chain consists of two components A and B, with Kuhn-lcngths 
I A and Ib and number of segments Na and Nb, respectively. There is no excluded- volume interaction and the solution 
is described by the linearized Navier-Stokes equation. The hydrodynamics is coupled to the chain dynamics by the 
no-slip boundary condition. 

Extending to our case the cluster expansion theory§0H§@ previously used to study the hydrodynamics of suspensions 
of polymers and spheres, we developed explicitly the necessary theoretical elements and we obtained analytically and 
calculated numerically the following physical quantities of interest: the diffusion coefficient Dab-, the friction coefficient 
Cab and the intrinsic viscosity [t/ab]- The results have been displayed graphically as functions of the dimensionless 
variables characterizing the problem, x = -jf- and t = yj. 

Even for the simple system studied here, an interesting behavior of Dab and implicitly of the other quantities, 
emerged. The main finding of the present work is that by adjusting x (the molecular mass ratio) and t parameters we 
can control the process of separation of diblock copolymers (or, more generally, combined objects) under hydrodynamic 
flow. To illustrate this aspect we constructed separation curves by taking the difference A(£>ab) = ^ab(^j^i) — 
Dab^j^) and representing it as a function of x and the pair t\ — ti, which allows us to identify both the absolute 
maxima and the ranges of x that maximize the difference A(£>ab) for given t\ and £2- As an immediate extension of 
the present work, an application of our approach to the process of electrophoretic separation is in progress. 

The formalism presented is applicable in more general contexts not limited to Gaussian chains. For example, 
when excluded-volume or Coulombic interactions are present we can capture their effects — in certain conditions — by 
renormalizing the Kuhn- lengths I a and Ib to some effective values. As an additional observation, the usual Kirkwood- 
Riseman approximation is not necessary. But, rather surprisingly, our calculations have shown rigorously that one 
can still safely apply it. 
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APPENDIX A: COMPUTATION OF (exp[ik(S s - S s /)]) 



Consider.-a, Gaussian polymer chain consisting of N links each of mean-square length Ik, with k <E l,N. Following 
YamakawaJlJ we can construct the probability distribution function P(Si — Sj) of the distance between two beads i 



and j located at and Sj , respectively, about the center-of-mass of the chain as 



2tt(Z 2 



3/2 



exp 



3IS, 



s,l 2 



2(1*) 



(Al) 
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with (l 2 ) — Ylm>j+i Then, (exp[ik ■ (S, — Sj)]) is the generating function of the moments of the P(Si — Sj) 
distribution and is given by 



M[k] = (exp[ik • {Si - Sj)}) = cxp(--k 2 (I 2 )). 



(A2) 



In our problem we have only two types of links, A and B, with average lengths I a and Ib and the calculation of the 
structure factor averaged over the distribution of beads yields: 



k 2 \s-s'\l A ,B 



if s, s' < La or s, s' > La; 



e Oc.(S s -S s ,)) = I e -lk 2 [(L A -s)l A + {s'-L A )l B ] jf g < La and s / > La . 

e -i k 2 [(L A -s')l A + {s-L A )lB] if g > La and s / < La . 



(A3) 



In the expressions above we introduced the arclength variables s and s' instead of the numbers i, j (e.g. s = Ha if 
s < La) and La = NaIa; Lb = NbIb- 



APPENDIX B: COMPUTATION OF THE MATRIX ELEMENTS G D 



As previously shown, the Fourier coefficients of the Oseen tensor are given by 



G, 



dk 1 - kk 



(2tt)3 r, k 2 JJ 



ds ds' 



(Bl) 



x (exp[ik ■ (S s - S s /)]) exp[-^ps + ^p's'] 



Introducing the structure factor from (A3) and performing the angular integral utilizing Jdflk (1 — kk) = we 
arrive at 



1 



W o 1 



dk- 



ds ds' e -^ 2 \s- s '\i A 







dsds'e-* k |s ~ s l' fl 

+ [ LA [ L dsds' e -\k 2 [(L A -s)l A + {s'-L A )l B ] 

Jo Jl a 

+ [ L dsds' e~* ^[(L A -s')l A + { s-L A )l B ] 

Jl a Jo 



x exp 



2z7r 1m . . 

--PS+—PS 



a. 



(B2) 



The evaluation of the integrals goes as follows. First we compute the Gqo factor: 



Gqo — 



16 1 

■^V^Vo IaIb 



R :5 



3 / l B 



>3 / l A 



-l)+Rf ? (-r i -l) + (R-a + R s ; 



.3/2 



1, 



(B3) 



where Ka,b — {La,b^a,b /G) 1 ^ 2 are the radii of gyration of the A and B components of the copolymer and the integral 
J™dx [£- £(1 - e-' 2 )] - was US ed. 

Next, making p = p' =/= and integrating over s and s' yields 



G PP = q o / dfc < 
37r 2 ?7o J 



2L\ (kR A ) 2 



2{L-L A f{kR B f 



( fcRs )4 + ^ 
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L A (L-L A ) 



x [1 — exp(- 



2inL 



^±p-k 2 R%)] 
L A {L-L A ) 



\2 I 2ttL a 



p][(kR 



i2 I 2tt(L — La) „ 



[l-exp(— ^p-fc 2 Ri)] 



;[l-ex P (-^p-fc 2 Ri)] 



(B4) 



x [l-expC^p-^Rl)]^. 



Integrating over in the complex plane the third and fourth factors of the integral vanish and we find that G pp is 
given by 



Gpp — L 



1 



7rV67r?7o y/L\p\ 



/IT L jT } 



1. 



(B5) 



The computation of the remaining elements G pp < | P ^ P ' and Go P is accomplished in a similar ma nne r and involves 
a rather lengthy calculus. A simplification occurs if noticing that the G pp ' operator given by (Bl) is hermitian 
(Gyp/ = G*, p ). Then we only need to calculate 



G, 



dk 



L A L 



(kR A ) 2 



3tt 2 ?7o Jo I ™{p' - P) (kR A ) A + ( ^ 



; 2 "j, A ip'-v) _ I 



L(L-L A ) 
i-k(j>' — p) 



1 - e ~ 



r — Hp p) 



(fcR B ) 4 + (M^kALp 



(B6) 



which yields the final result 



G PP ' — L 



2ir 2 Vfar]o y^L\p~\ i(p' ~ p) \ Vh ^/Ta 
Letting p' — > gives G p o and from hermiticity we obtain Gq p = G* . 



1 - e L 



2iz r lA (p'-p) 



(B7) 



APPENDIX C: EVALUATION OF THE AVERAGE (S s S s /) 



We consider a Gaussian diblock copolymer A-B consisting in N A links of Kuhn-lcngth l A and Nb links of Kuhn- 
lcngth Is, with N A + Nb = N. Let be the bond vector (along the segment) of the fc-th link. Our goal is to 
construct the bi-variate probability distribution function P(Si, Sj).t2l This can be done by expressing the positions 
and Sj of the «-th and j-th beads about the center-of-mass as 



N N 



Si = ^ik?k\ Sj = Y ^jkr k , (CI) 



k=l k=l 



where 



* tt = 6(i -k) + - 1; ¥ a - 9(j -k) + - 1 (C2) 



and Q(x) is the step-function. Also, we introduce the symmetric matrix 



Ca Cij 



(C3) 
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with its components given by 



1 / NA N \ 

7i2\ ( l A ^rnk^nk + l B y^nk I 

^ ' V k=l k=N A +l ) 



m,ne {i,j} 



(C4) 



and (; 2 ) = jf[N A l 2 A + {N-N A )l 2 B ] = l 2 A (x + t 2 {1 - x)) . 

With these notations, we can calculate in a well documented fashion0 the bi-variate probability distribution function 
of the Gaussian diblock copolymer A-B as 



P(S i5 Si) 



2tt (Z 2 ) y (detC) 3 / 2 



exp[-CiS 4 • Sj - C 2 Si ■ Sj - C 3 Sj • Sj], 



(C5) 



where 



C = 



2 (I 2 ) det C 



^ 3Cij 3C,;, 



(/ 2 )detC ° 2(l 2 )detC 



(C6) 



The generating function for the moments is equal to: 

1 



M[k,k'] = exp 



-- (I 2 ) (C 4l k 2 + 2Cyk • k' + C n V 2 ) 



(C7) 



The sought average (Si • Sj) results immediately from: 



(Si • Sj) = lim z 2 V k V k ,M[k,k'] = (I 2 ) Qj 

k— >0 



(C8) 



It is clear that the value of the CV,- factor depends on where the beads i and j are located with respect to the two 
sections of the copolymer. We have three cases: both beads in block A (C^j), both beads in block B (C^) or the 



i-th bead in A and the j-th bead in B (C^ B = C BA ). Inserting the values of ^(i,j)k fr° m E q. (|C2|) in Eq. (JC^) and 
carrying out the summation keeping only terms of order N, we obtain 



C~ = 



C» = 



;2 
M 



2A 



(Z 2 

iV 



2A 



+ t 2 



,iV 



;2 
M 



(Z 2 ) 

(/ 2 ) 



3 / 3 

[< 2 + (i-< 2 )x 3 ] + i 2 



1 — 3x 1 — x 



2N 



i 2 1 / 2 , 
2-V 2 A "'(' 3* 1 



r 

2N 
2 



2N 3 + 3 



iV A 3 



1 + ^x 2 - 



(C9a) 
(C9b) 
(C9c) 
(C9d) 



where we remind that x = ^ and t = j 3 -. Returning to the arclength variables s, s' (e.g. s = ii^ when s < La or 
s = NaIa + (i — Na)Ib for s > Na), we complete the computation of (S s • S S ') by inserting the CV,- factors in (|C 
In the three mentioned cases of distribution of the beads s (former i) and s' (former j) in the blocks A and B we find 



{Ss • S 5 



(S x • Sv) 



[s 2 + s' 2 - 2sNl A 0(s - s') - 2s'Nl A {\ - 6»(s - s')) 



2N 

+ 2{Nl A f [x{\ - t 2 )(l -x+ ^x 2 ) + U 2 

1 

2A 



(ClOa) 



[s 2 + s' 2 - 2Nl A [x{\ -t)- t9(s - s')]s 



-2Nl A [x{\ -t)+ i(l - 8(s - s'))]s' 
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/I +2 

+ 2{Nl A f ( -x 2 (l - t 2 ) + - + xt(l -t) + x 2 (l - t)' 
1 



(S s • S S ,) AB = — [s 2 + s' 2 - 2Nl A (t + x{\ - t))s' 

+ 2{Nl A ) 2 Qx 2 (l - t) 2 + xt(l -t) + j 



-x 2 (l-< 2 )(l--x) 



(ClOb) 



(ClOc) 



with (S, ■ S S ,) BA = (S, ■ S,') AB 



Having calculated these averages, we introduce them back in 



F / — F /H 

pp — pp 



(S s ■ S S ') exp[- 



2m 2m 



L 



■ps 



L 



p's'] t. 



and after evaluating the resulting integrals we arrive at the Fou rier c oefficients F pp /. We show the /oo, 
elements (where / pp < is the dimensionless part of F pp < from Eq. ( 3.31 )): 



(Cll) 
fpp and f 0p 



foo(x,t) 
fpp(x : t) 



fo P (x,t) 



(l-tf(l-x) 2 x 2 (t 2 + x-t 2 x) 



4p 3 7T 3 



-2i P 7TX 

e'+'-'' (— t — x + tx 



> 2 (- 



■i t + i e t +"- tx t + i t —ie t + i - tx t 



_l_4 e t+x-tx put — ix + ie*+*-** x + 4e'+'-i' pirx 

+ 2 i t x — 2 i e tx — it 2 x + i e f +'"- ti: t 2 x — 4 e put 2 x 



1 



(-1 + t) (-t-x + ix) (it 2 -ie'+"-'* i 2 



4tt 3 p 3 

+2 i < x — 2 i e*+ a =- ta: i x — 2 e'+"-" it ptx — 2it 2 x + 2i e T +^ T ^ t 2 x 
+2 e '+*-'■ 7r p i 2 a; + 2 1 7r 2 p 2 t 2 cc + i x 2 - i e'+*-'' x 2 

2 \ tt p x q Q O O 

—2 e*+ x - tx 7rp x +2i7r p x 

-2 i < x 2 + 2 i e t+J-tl tx 2 +Q e tt ptx 2 + it 2 x 2 - i e '+'*-'■ i 2 x 2 

_4 e •+*-'■ 7T p i 2 x 2 - 4 i 7T 2 p 2 t 2 x 2 + 2 e '+*-'" 7r p x 3 - 2 i 7r 2 p 2 x 3 
_ 4 e t+x-tx tt ptx + 2e t +*- t * irpt x r +2iir p t x 6 



(C12a) 
(C12b) 



(C12c) 



Unfortunately, f vv >( x, t) i s too long to be reproduced here. Limits of / pp '(x,i) when x = {0, 1} and £ = 1 were quoted 
in the text (see Eq. (3.32)). 
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Diffusion Coefficient of the A-B chain 

N A = # of segments in A; N = total # of segments in A-B; t = 1 B /1 A 

10.0 r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r ~ 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



x = N A /N 

FIG. 1. Diffusion coefficient Dab/D as a function of x = Na/N and t = Ib/Ia (D is the diffusion coefficient of a 
homogeneous chain with length Ljv = NIa) '■ solid line t = 1/10; dotted line t = 1/8; long dashed line t = 1/6; dashed line 
t = 1/4; dot-dashed line t = 1/2; thick solid line t = 1; o-long dashed line-o t = 2; o-dashed line-o t = 4; o-solid line-o t = 10 



Diffusion Coefficient of the A-B chain 

N A = # of segments in A; N = total # of segments in A-B; t = 1 B /1 A 




x A = N A /N 

FIG. 2. Diffusion coefficient Dab/Da as a function of x = Na/N and t = Ib/Ia (Da is the diffusion coefficient of a 
homogeneous chain with length La = NaIa): solid line t = 1/10; dotted line t = 1/8; dot-dashed line t = 1/6; dashed line 
t — 1/4; long dashed line t — 1/2; thick solid line t = 1; o-long dashed line-o t = 2; o-dashed line-o f = 4; o-dot-dashed 
line-o £ = 6; o-dotted line-o £ = 8; o-solid line-o t = 10 . 
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Separation Curves 



N A = # of segments in A; N = total # of segments in A-B; t = 1 B /1 A 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

x = N A /N 



FIG. 3. Separation curves as A(D A b/D a )(x) = ^f-{x,ti) - ^f(x,t 2 ): solid line [ti = 1/8; t 2 = 1/6]; dotted line 
[ti = 1/6; t 2 = 1/4]; long dashed line [ti = 1/4; t 2 = 1/2]; dashed line [ti = 1/2; t 2 = 2]; dot-dashed line [ii = 2; t 2 = 4] . 

Friction Coefficient for the A+B chain: t<l 



N A = # of segments with 1 A ; N = total # of segments in A-B; t = 1 B /1 A 




' 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



x = N A /N 



FIG. 4. Friction coefficient Cab/C as a function of x and t (£ is the Kirkwood-Riseman (K-R) friction coefficient for 
a chain with length Ln = NIa)- Thin lines show the results obtained when using only the diagonal elements of G pp i 
for the t values corresponding to the symbols: thick solid line t = 1; o t = 1/2; □ t = 1/4; A t = 1/6; \/ t — 1/10; - 
Diagonal G pp i. 
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Friction Coefficient of the A-B chain: t>l 



N A = # segments in A; N = total # of segments in A-B; t = 1 B /1 A 




' 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



x = N A /N 

FIG. 5. Friction coefficient Cab/C as a function of x and t (£ is the K-R friction coefficient for a chain with length 
Ljv = NIa )• Thin lines show the results obtained when using only the diagonal elements of G pp / for the t values 
corresponding to the symbols: thick solid line t = 1; o t = 2; □ t = 4; A t = 6; O t = 8; V * = 10; — Diagonal G pp i . 

Change in viscosity due to the A-B copolymer 

N A = # of segments in A; N = total # of segments in A-B; t = 1 B /1 A 




x = N A /N 



FIG. 6. Change in the viscosity St/ab/Sti as a function of x and t (Sr/ is the K-R non-free draining result for the change 
in viscosity due to a polymer with length L = NIa)- Thin lines show the results obtained when using only the diagonal 
elements of G pp i for the t values corresponding to the symbols: thick solid line t = 1; o t = 2; □ t = 3; A t = 4; O < = 5; 
V t = 6; — Diagonal G pp t . 
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Change in viscosity due to the A-B copolymer (K-R approx.) 

N A = # of segments in A; N = total # of segments in A-B; t = l B fl A 




" 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

x = N A /N 

FIG. 7. Change in the viscosity SrjAs/SrjA as a function of x and t using the Kirkwood-Riseman approximation (Sr/A 
is the K-R non-free draining result for the change in viscosity due to a polymer with length La = NaIa)'- thick solid line 
t = 1; dotted line t — 4/5; long dashed line t = 3/5; dot-dashed line t — 2/5; dashed line t — 1/5; solid line t — 1/10 . 
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